Turbulent front speed in the Fisher equation: dependence on Damkohler number 
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Direct numerical simulations and mean-field theory are used to model reactive front propagation 
in a turbulent medium. In the mean-field approach, memory effects of turbulent diffusion are taken 
into account to estimate the front speed in cases when the Damkohler number is large. This effect 
is found to saturate the front speed to values comparable with the speed of the turbulent motions. 
By comparing with direct numerical simulations, it is found that the effective correlation time is 
much shorter than for non-reacting flows. The nonlinearity of the reaction term is found to make 
the front speed slightly faster. 



I. INTRODUCTION 

It is well known that the propagation speed of a flame 
front is greatly enhanced if a mixture of fuel and oxy- 
gen is in a turbulent state. This topic of turbulent pre- 
mixed combustion was pioneered by Damkohler |l| some 
70 years ago and is reviewed extensively in recent liter- 
ature 2-4]. In spite of its importance, the question of 
burning velocities in a turbulent medium continues to be 
of major importance even today 

Much of the current work is based on the original 
Damkohler paradigm for premixed combustion. He dis- 
tinguishes two regimes, namely those of large-scale and 
small-scale turbulence. For the purpose of the present 
paper it is useful to base this distinction on a compari- 
son of the mean turbulent flame width with the scale of 
the energy-carrying eddies @j. In the small-scale turbu- 
lence regime, also referred to as the distributed reaction 
zone regime, the turbulent flame speed is computed using 
a formula where the microscopic diffusivity is replaced by 
the sum of microscopic and turbulent diffusivities. This 
is possible because there is good scale separation. This 
implies that the turbulent front thickness (i.e. the thick- 
ness of the flame brush) is much broader than the scale 
of the turbulent eddies. This regime is characterized 
by small Damkohler numbers. In the opposite case of 
large Damkohler numbers, the turbulent front thickness 
is smaller than the scale of the turbulent eddies and can 
therefore no longer be described by turbulent diffusion. 
This regime is characterized as that of large-scale turbu- 
lence. In this case the turbulent front speed reaches its 
maximal value that is given by the rms velocity of the 
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turbulence in the direction of front propagation. 

The regime of large-scale turbulence is subdivided fur- 
ther into regimes of corrugated and wrinkled flamelets, 
depending essentially on the ratio of Reynolds number to 
Damkohler number, which is also related to the Karlovitz 
number. When the Reynolds number is small compared 
with the Damkohler number (small Karlovitz number), 
the flame front is merely wrinkled, but for large Reynolds 
numbers (large Karlovitz number) it becomes corrugated 
and can consist of isolated flamelets detached from other 
parts of the front. In the present paper we will mainly be 
concerned with the flame speed rather than the question 
of whether the flame front is wrinkled or corrugated. 

In turbulent combustion, the averaged flame speed, St, 
is usually normalized by the corresponding laminar flame 
speed, sl, and one is interested in the dependence on 
the normalized turbulent velocity, v'. For the regime of 
large-scale turbulence, the speed-up ratio of turbulent to 
laminar flame speed is given by the geometric increase of 
the wrinkled surface area of the flame front. Damkohler 
assumed that the increase in surface area is proportional 
to the ratio of the turbulent velocity of the eddies to the 
laminar flame speed. This leads to the expectation that 
the dependence of st on v' is given by 



st/sl = 1 + v'/s L . 



(1) 



This equations captures the expected limiting cases that 
st should not become larger than v' and that st = sl 
in the absence of turbulence, i.e., for v' = 0. However, 
unsatisfactory agreement with measurements motivated 
the search for other dependencies. For example, Pocheau 
Q derives the more general formula 



st/sl = [1 + (v'/s L ) n ] 



l/r 



(2) 
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where n is a parameter. This formula obeys the afore- 
mentioned limiting case for any value of n. Pocheau 
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8] contrasts the formula with another one proposed by 
Yakhot 0, 



st/sl = exp (v'/s T ) 



(3) 



where st < v' 
given by 



for v' — > oo. Yet another fit formula is 



s T /s L = 1 + C w (v'/slY 



(4) 



with fit parameters Cw and m = 0.7 10]. Both (j3|) and 
(|4]) have a front speed less than v' for v' oo, provided 
to < 1 in Eq. Q. As can be seen from Fig.[TJ the different 
proposals for the front speed are quite similar, making 
it difficult to use measurements to distinguish between 
them. Furthermore, realistic descriptions of flame prop- 
erties are hampered by the fact that feedback on the flow 
by the actual combustion process depends to the specific 
case and is not easy to model. The feedback on the flow is 
therefore usually ignored. It might therefore be useful to 
return to a simple model of front propagation that can 
be treated in more detail and to address the unsettled 
question regarding the different proposals in Eqs (J)-® 
for the dependence of St on v' . Following Kerstein [ll| , 
we consider here the Fisher equation, which is also known 
as the Kolmogorov-Petrovskii-Piskunov (KPP) equation 
[l9j . An important difference to earlier work is the fact 
that we solve this equation in the three-dimensional case 
in the presence of a turbulent velocity field. 

The Fisher or KPP equation is a simple scalar equation 
that possesses propagating front solutions. This equa- 
tion is familiar in biomathematics [l2j as a simple model 
for the spreading of diseases. It has also been amended 
by an advection term to describe the interaction with 
a turbulent velocity field in one [li| and multiple [l5[ 
dimensions, the effects of cellular flows and the scal- 
ing of the front thickness [TtJ • Furthermore the equation 
has also been modified to account for different interact- 
ing species, that can be used to model the spreading of 
auto-catalytically polymerizing left and right handed nu- 
cleotides [18[ . Given that C is a passive (albeit reacting) 
scalar, the Fisher equation does ignore any feedback on 
the flow and is therefore well suited to help clarifying 
questions regarding the relation between st and v' . 

In the present paper we consider both direct numerical 
simulations (DNS) of this equation in three dimensions as 
well as its averaged form where the effects of turbulence 
are being parameterized by a non-Fickian diffusion equa- 
tion. Such an equation allows for the ballistic spreading 
of a passive scalar concentration on short time scales, 
which is expected to be important when the front prop- 
agates at a speed comparable to that of the turbulence 
itself. 



II. THE FISHER EQUATION 

A simple model of front propagation is the Fisher equa- 
tion which, in the simplest case, is a one-dimensional 
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FIG. 1: Comparison of different expressions for the normal- 
ized front speed, st/sl, as a function of the turbulent velocity, 
u'/sl- The labels n = 1 and n — 2 refer to Eqs. |T]) and ([2]), 
while Y88 and W85 refer to Eqs. © and ©. 



partial differential equation [12|, I19H21 



dC _ C t , _C_ 



D 



, d 2 C 
dx 2 ' 



(5) 



for the concentration C. Here, r c is the chemical reaction 
time, D is the diffusivity, and Co is some saturation value 
above which further growth is quenched. Equation ([5]) 
corresponds to an autocatalytic reaction were a reactant 
TZ yields a product V at a rate A; that is itself proportional 
to the concentration of the products, [V], i.e. 



K 



V with k = [V]/t c C q . 



(6) 



This can then also be written as V+TZ —> IV. Saturation 
of the product concentration, C = [V], results from the 
fact that the total mass is conserved, i.e. [TZ] + [V] = 
Co = const. The evolution equation for the concentration 
C = [V] is then given by Eq. ©. 

This equation has two solutions, an unstable solution, 
C = 0, and a stable one, C — Cq. The diffusion term 
seeds the neighboring regions that are in an unstable 
state, which promotes the rapid transition from C = to 
C = Co. This leads to the propagation of the transition 
front in the direction down the gradient of C with a front 
speed [H} 



(7) 



where the subscript L refers to the laminar front speed. 

In many cases of practical interest the diffusion co- 
efficient D is rather small and is hardly relevant when 
there is rapid advection through fluid motions. In that 
case the governing equations become advection-reaction- 
diffusion equations. This can be written as 



at t c V Co 



L»V 2 C, (8) 
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where U is the flow speed. If the flow is turbulent and has 
zero mean, there can be circumstances where the average 
concentration C can be described by an equation similar 
to Eq. (|5|) , but with C being replaced by the mean value 
C, and D being replaced by some turbulent diffusivity 
D t , i.e. 

dC C ( C\ „ d 2 C 

where Dt = D + D t is the total (i.e. the sum of micro- 
scopic and turbulent) diffusivity. We have here assumed 
that the mean concentration shows a systematic varia- 
tion in the x direction and have thus assumed averaging 
over the y and z directions, so C{x,t) can be described 
by a one-dimensional evolution equation. 

Given the similarity between Eqs. (O and ©, one 
would expect that in the turbulent case with appropri- 
ate initial conditions the effective turbulent propagation 
speed st of the front can still be described by an expres- 
sion similar to Eq. ([7]), but with D being replaced by -Dt, 
i.e. st = 2 yj Dt / r c . A useful estimate for the turbulent 
diffusivity is D t — u vms /3k t , where k t is the wavenumbcr 
of the energy-carr ying eddies and u rms is the rms velocity 
of the turbulence [13]. Thus, for D t ^> D, the effective 
value of st is expected to be 2(u rms /3r c fc t ) 1/ ' 2 - On the 
other hand, one cannot expect the front speed to increase 
indefinitely with decreasing r c . Indeed, one would not ex- 
pect st to exceed the rms velocity of the turbulence in 
the direction of front propagation. Following common 
practice, we denote it by v' . Under the assumption of 
isotropy, v' is related to the three-dimensional rms veloc- 
ity by v' = u rms /V3. 

An important nondimensional measure of r c is the 
Damkohler number, which is the ratio of the turnover 
time, (u lms kt)~ , to r c . This number is here defined as 

Da= (rctwsfct)- 1 . (10) 

Note that our definition of Da is based on the wavenum- 
ber kt rather than the scale 2w/kt, which would have 
reduced the numerical value of Da by a factor of 2ir. For 
small values of Da we expect st ~ 2t/ Da 1/2 , while for 
large values one expects st « v' 0. Thus, a more general 
formula is expected to be 

s 2 t = «' 2 /(Da), (11) 

where /(Da) increases linearly with Da for Da <C 1 and 
/(Da) w 1 for Da 3> 1. This saturation behavior can also 
be interpreted as a reduction of the effective value of r c 
[24j . An important goal of this paper is to determine the 
form of the function /(Da). 

III. NON-FICKIAN DIFFUSION 

The Fickian diffusion approximation made in Eq. @ 
for the mean concentration C becomes invalid if C varies 



rapidly in time, and in principle also in space. This is 
indeed expected to be the case when Da 3> 1. For rapid 
time variations, Eq. ([9]) attains then an extra time deriva- 
tive and takes the form [2o| 

d 2 C dC V / C\ „ d 2 C 

which is a damped wave equation with relaxation time 
r and an additional reaction term. The presence of the 
nonlincarity in the reaction term leads to an additional 
contribution in the C equation which has here been ig- 
nored (see Appendix lAl for a more consistent treatment). 

Without the reaction term, Eq. (|12j) is also known as 
the telegraph equation. This equation emerges naturally 
when computing turbulent transport coefficients using 
the r approximation (23[. Evidence for the existence of 
the wave term has been found from isotropic forced tur- 
bulence simulations [25]. A non-dimensional measure of 
t is given by the Strouhal number, 

St = Tu lms k t = Tu 2 ms /3A, (13) 

where the first equality is useful for turbulence simula- 
tions where ru rms /c t is readily evaluated, while the second 
equality is useful for the mean-field model, where k t does 
not appear explicitly and D t and u rms are given. 

Using DNS of forced turbulence with a passive scalar, 
the value of St has been determined to be around 3 by 
relating triple corrections to quadratic ones [25|. Al- 
though we consider the value of St as being fairly well 
constrained, we do consider below a range of different 
values. 

The purpose of this section is to study solutions of 
Eq. (|12p that can then be compared with DNS of the 
Fisher equation coupled with the Navier-Stokes equations 
for obtaining a turbulent velocity that enters Eq. JS} . We 
consider first the case where D is negligible and solve 
Eq. (TT2")) for different values of Da in a one-dimensional 
domain that was chosen long enough so that the front 
speed can be determined accurately enough. We use a 
numerical scheme that is second order in space and third 
order in time [26j]. In some cases a resolution of 2 15 
3 x 10 4 mesh points was necessary. 

We study first the dependence of the front speed on Da 
for a range of different values of St and D <C D t . Here, 
st is determined by differentiating the concentration in- 
tegrated over the whole domain, 

and approximating the asymptotic front speed with the 
value at the time when the front has reached the other 
end of the domain. This quantity is also known as the 
reaction speed. This is indicated by C reaching a small 
fraction (e.g. 1CU 6 ) of Co- The result is shown in Fig. [2] 
For small values of Da, the front speed is independent of 
the value of St and we reproduce the anticipated result, 
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FIG. 2: Dependence of the front speed of solutions of Eq. (|12[) 
on Da for different values of St and Pe = oo. The lines repre- 
sent fits given by Eq. (|15[) . 



FIG. 3: Dependence of the front speed of solutions of Eq. (|12[) 
on Da for different values of Pe and St = 3. The lines repre- 
sent fits given by Eq. (fl8|) . 



i.e. /(Da) = 4 Da. For larger values of Da the front speed 
reaches eventually a constant value. However, the limit- 
ing value depends on St. Our results are well reproduced 
by the fit formula 



= /(Da, St) 

o J* 



4 Da 



1 + 3 St Da 



(15) 



This formula obeys the anticipated limiting behaviors for 
small and large values of Da, provided St ss 4/3. The nu- 
merically determined data agree quite well with Eq. (|15p . 
However, in some cases the numerical data are somewhat 
uncertain and depend also slightly on resolution and do- 
main size. 

In the DNS presented below, where the numerical res- 
olution is still limited, the value of D is often not negli- 
gible. Its value is characterized by the Peclet number, 



Pe = u xma /Dk t « 3D t /D. 



(16) 



For small values of Da, the expr ession for the front speed 
should be st = 2 \J (D t + D)/t c , which can then be writ- 
ten as 



s T /v' = 2^(1 + 3P6- 1 ) Da (forDa<l). (17) 

For larger values of Da we solve Eq. (fT2"j) numerically. A 
good fit formula for / is given by 



/(Da, St, Pe) « 4Da 



3 

Pe 



1 



1 + 3StDa 



(18) 



In Fig. [3]we compare the fit formula with the numerically 
obtained front speeds for different values of Pe, keeping 
St = 3. The agreement is again quite good. 

In turbulent combustion it is customary to plot the 
normalized front speed, st/sl 7 as a function of the nor- 
malized turbulent velocity, u'/sl- Using our definitions 



of Da and Pe in Eqs. (ITU1) and (|TB)) . respectively, we have 
v'/sl = (Pe/12Da) 1 / 2 and find 



3St/4+eu'/s L 



1 1/2 



(19) 



where we have defined e = fct^F with £f = (r c D/12) 1 / 2 
being a measure for the laminar flame thickness. Note 
that e can also be expressed in terms of Da and Pe via 



(12DaPe) 



-1/2 



(20) 



A more familiar quantity is the ratio £/£f — e , where 
£ = k^ 1 is the typical eddy scale. Even if we can assume 
the value of St to be given, e is not a fixed quantity. It 
is therefore clear that there cannot be a unique relation- 
ship between st/sl and v'/s^. Instead, there must be a 
family of solutions depending on the value of e; see Fig. 2) 



IV. DNS OF THE FISHER EQUATION 

We now consider DNS of Eq. © where U is obtained 
by solving the Navier-Stokes equation for an isothermal 
gas with a forcing term that is S correlated in time. The 
forcing function consists of plane waves whose wave vec- 
tor is random and its length is within a narrow window 
around some mean forcing wavenumber kt- Since the gas 
is compressible and the density p is not constant, Eq. flSJ) 
now takes the form 



dC_ 

dt 



u -vc = 



c 



1 - 



c_ 



which we solve together with the momentum and conti- 
nuity equations, 



8U 
~dt 



-U-VU-c 2 s Vlnp + f + F lolc 



(22) 
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TABLE I: Summary of the runs discussed in this paper. 



Run A 


He 


1 1 _ 

rc 


Da 


Ka fct/fci 


St/V 


st/sl 


v'/sl £/£f 


A 1 


1 


117 


117 


0.2 


1 A £. 
14.0 


O.U 


0.97 


7.50 


"7 7£ 
1.10 


10. U 


A O 
Az 


1 


115 


115 


0.5 


4.0 


. 1 


1.60 


7.03 


A QQ 


of; q 
ZD.o 


A • >. 
A.i 


i 

i 


122 


49 


1.6 


O.O 


□ . 1 


2.33 


3.77 


1.01 


oU.4 


A4 


1 


121 


12 


4.8 


4.8 


5.1 


3.55 


1.63 


0.46 


26.3 


A5 


1 


120 


4 


15.9 


3.6 


5.1 


7.29 


1.16 


0.16 


30.3 


Bl 


4 


38 


513 


1.6 


3.6 


1.6 


1.88 


9.76 


5.19 


98.9 


B2 


2 


41 


164 


4.9 


3.6 


1.6 


2.21 


3.69 


1.67 


98.8 


B3 


2 


165 


165 


4.9 


3.6 


1.6 


2.52 


4.22 


1.67 


98.9 


B4 


2 


41 


16 


4.9 


36.4 


1.6 


2.80 


1.47 


0.53 


31.3 


B5 


2 


40 


16 


50.3 


3.6 


1.6 


7.25 


1.19 


0.16 


98.9 



FIG. 4: Dependence of st/sl and v' /sl for e = (solid line), 
0.3 (dotted), 1 (dashed), and 3 (dash-dotted). Note that there 
is no unique relationship between st/sl and h'/sl- 

| = -V.pt/, (23) 

where F forcc = zy(V 2 t/ + ±VV • (7 + 2SVlnp) is the 
viscous force, S = \ [VU + (VC/) T ] - |lV ■ U is the 
traceless rate of strain tensor, I is the unit matrix, v is 
the molecular diffusivity, and c s = const is the isothermal 
sound speed. The forcing function / is of the form 

f(x, t) = Re{Nf k{t) exp[i*(t) • x + ty(t)]}, (24) 

where x is the position vector. The wave vector k(t) and 
the random phase — it < <fi(t) < tt change at every time 
step. For the time-integrated forcing function to be inde- 
pendent of the length of the time step 5t, the normaliza- 
tion factor N has to be proportional to 5t~ x / 2 . On dimen- 
sional grounds it is chosen to be N — foc s (k t c s / 'St) 1 / 2 , 
where f is a nondimensional forcing amplitude. The 
value of the coefficient /o is chosen such that the maxi- 
mum Mach number stays below about 0.5. Here we 
choose /o = 0.02. We force the system with nonhelical 
transversal waves, 

f k = (kx e) /Vfc 2 -(fc-e) 2 , (25) 

where e is an arbitrary unit vector not aligned with k; 
note that \ fk\ 2 = 1- 

In the x direction we use periodic boundary conditions 
for U and p and dC/dx — for C, while we use peri- 
odic boundary conditions in the y and z directions. The 
simulations were performed with the Pencil Code [22j], 
which uses sixth-order explicit finite differences in space 
and a third-order accurate time stepping method 26] . 

For the calculations we use units where k\ — c s — 
po = 1. However, most of the results are presented in 
an explicitly non-dimensional form by normalizing with 
respect to relevant quantities such as the rms velocity of 
the turbulence or the turnover time. Our simulations are 
characterized by several non-dimensional parameters. In 



addition to the values of Da and Pe, defined in Eqs. © 
and p6p . respectively, there is the Schmidt number, Sc = 
v/D. In those cases where the Damkbhler number was 
large, we had to increase the value of D in order to resolve 
the flame front. This was done by decreasing Sc to values 
below unity. The degree of scale separation is given by 
the ratio kt/k\. 

V. RESULTS 

In the following we present results for uniform aspect 
ratio, A — L x /L y = 1 with k t /ki — 5 (series A), and 
A = 2 or 4 with kt/k\ = 1.6 (series B). Our runs of 
series A and B are summarized in Table HI The resolu- 
tion in the y and z directions is always 256 2 meshpoints, 
but it is larger in the x direction in runs where the as- 
pect ratio A is larger than unity. In Fig. [5] we show the 
concentration C on the periphery of the computational 
domain at different times for r c = 3/c s ki, which corre- 
sponds to Da = 0.5; see Table HI One sees clearly how 
the front spreads and propagates in the negative x di- 
rection. The front speed is determined in the same way 
as for the mean-field model, i.e. using Eq. (TH| . except 
that C is computed from the actual C. This can also be 
formulated as a volume integral, 

^ ^ L x L y dt J Co ^ ^ 

In Fig. |5] we show examples of the evolution of the mean 
concentration and the instantaneous front speed as func- 
tions of time for series A. The resulting ratios st/v', 
s t/sl, ^'/ s l, and £/£f are summarized in Table U for 
series A and B. 

In most of the cases considered in this paper, the value 
of Pe is not in the asymptotic regime. It might therefore 
be sensible to compare the relative front speed, st/v' 
against the function /(Da, St, Pe). This is done in Fig. [71 
where we show the non-dimensional front speed, St/v', 
versus /(Da, St, Pe), for three values of St using values 
of Da and Pe, as evaluated from Eqs. ([TUf and (|16l) . Sur- 
prisingly, the best fit is obtained for rather small values 
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FIG. 6: Mean concentration and the instantaneous front 
speed as functions of time. 



of St of 0.03. This suggests that, for the present applica- 
tions, the relevant value of r is much smaller than in the 
case of a non-reacting passive scalar. 

Next, we plot st/V versus Da for different values of 
Pe; see Fig. fusing the previously inferred value St=0.03. 
The data points from the DNS tend to lie between the 
curves for Pe = 1 and 10, even though most of the actual 
values of Pe are beyond Pe=10. This too suggests some 
inconsistency between the DNS and the mean-field de- 
scription in terms of the telegraph equation. Finally we 
plot the DNS results in a state diagram of St/sl versus 
v'/sl using St=0.03; see Fig. |H1 The data lie between 
the theoretical curves for £/£p = 10 and 100, which is 
roughly in agreement with the values given in Table HI 
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FIG. 7: Relative front speed as a function of / for three 
values of St. The squares indicate runs where the fluid is 
at rest and the front is moving through the domain while 
the asterisk denote runs with an inlet velocity chosen such 
that the front is approximately stationary within the do- 
main. The solid line gives the theoretically expected result, 
st/v' = /(Da, St, Pe) 1 / 2 . Note that the best agreement with 
the theoretical values is achieved for St=0.03. 



VI. CONCLUSIONS 

In the present work, the Fisher equation has served as 
a simple model equation for front propagation in a tur- 
bulent flow. The model has similarities with turbulent 
combustion, but is much simpler. Nevertheless, it is clear 
that even this simple model harbors surprises that one 
might have overlooked under more complex conditions. 



FIG. 5: Visualization of the concentration C on the periphery of the box at different times for Run Al. Here, T — 
is the turnover time. 
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a non-Fickian diffusion law for the turbulent diffusion 
and to retain the time derivative in the expression for 
the concentration flux. Earlier work did already confirm 
the principal validity of this approach and resulted in an 
estimate for the relevant relaxation time, which is charac- 
terized by the Strouhal number. The current work shows 
that the best fit to the simulation data can be achieved 
with a Strouhal number that is as small as 0.03, which is 
about 100 times smaller than the earlier determined value 
for passive scalar diffusion in forced turbulence. This dif- 
ference is connected with the presence of a reaction term 
in the evolution equation for the passive scalar concen- 
tration. 



Da 



FIG. 8: Relative turbulent front speed versus Da. The 
squares indicate runs where the fluid is at rest and the front 
is moving through the domain while the asterisk denote runs 
with an inlet velocity chosen such that the front is approx- 
imately stationary within the domain. For the latter, big 
asterisks denote cases where Pe > 10. The lines give the the- 
oretical expectations for St = 0.03 and Pe=l (solid line), 10 
(dotted), and 100 (dashed line). 
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FIG. 9: Turbulent front speed versus turbulence intensity 
for e = 0.1 (solid line), 1 (dotted), and 10 (dashed) using 
St = 0.03. The squares indicate runs where the fluid is at 
rest and the front is moving through the domain while the 
asterisk denote runs with an inlet velocity chosen such that 
the front is approximately stationary within the domain. For 
the latter, big asterisks denote cases where Pe > 10. The lines 
give the theoretical expectations. 



Using three-dimensional simulations we have been able 
to compare with the associated mean-field model. For 
small Damkohler numbers, the effective front speed can 
be approximated by replacing the diffusivity by a tur- 
bulent value. However, for Damkohler numbers larger 
than unity, this simple procedure fails, because it would 
suggest front speeds that exceed the characteristic speed 
of the turbulent eddies. A simple remedy is then to use 



Appendix A: Mean-field effect of the reaction term 

In order to assess the effect of neglecting the reaction 
term in the analysis presented above, we present now 
a simple mean-field theory for the Fisher equation using 
the r equation. We start with the passive scalar equation 
with a reaction term as given by Eq. ©, split C = C + c 
and U = U + u into mean and fluctuating parts, neglect 
the molecular diffusion term for simplicity, and define the 
mean concentration flux !F = uc and the mean squared 
concentration, % — c 2 , so the equation for the mean 
concentration is 



1 - 



C_ 

C 



H 

t c Cq 



, (Al) 



so the equation for the fluctuations is, to linear order in 
the fluctuations, 



9C _ ,== T=7x c 

= -V-(C/c + wC) + - 

Ot T c 



2C 
Co~ 



(A2) 



where the dots denote higher order terms for which we 
shall adopt a general closure assumption. Next, we derive 
evolution equations for J- and H, ignore a mean flow for 
simplicity, and assume V • u = 0, so we have 



m 

dt 



-2J-VC + 2 



u 



2C 
Co 



H 

T 



(A3) 
(A4) 



In Eqs. (|A3J) and (|A4j) we can write the last two terms 
as —T/tjt and —%/t-h, respectively, where 



1 



i 

tuQS) 



1 




(- 


r 


Tc 




1 


2 




r 







2C 

2C 
Co 



(A5) 
(A6) 



On sufficiently long time scales we may ignore the time 
derivatives in Eqs. (IA3|) and (IA4[) . so we arrive at closed 
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TABLE II: Dependence of st/v' without and with H in a 
model for Pe = 10. Note the slight increase of st/v' when H 
compared to the case where it is neglected. 



Da 


st/v' (without H) 


s T /V (with H) 


0.10 


0.25 


0.25 


0.30 


0.44 


0.47 


0.50 


0.59 


0.65 


0.61 


0.66 


0.73 



expressions for T and H, that we insert into Eq. (|Al|) . 
so we obtain 

f + «.*7=!(i-§) +Dl VC, (A7) 

where U C (C) is a new effective advection speed and Dt — 
D + D t is again the sum of turbulent and microscopic 
diffusivities with 

[/ C (C)=2A — D t (C)=rr(C)v' 2 . (A8) 

T~c Gq 

One may expect that the term U c slows down the prop- 
agation speed of the front, because it is directed up the 
concentration gradient. Note that the sign of the U c term 
is opposite to that of a similar term in the so-called G 
equation |2UlCj of turbulent front propagations, which is 
however not an equation for the flame brush, but for the 



detailed position of the wrinkled flame front (at G = 0) 
with an advection speed that is given by u ~ s^n, where 
n = VG/|VG| is a unit vector normal to the flame front, 
but it enters with a minus sign and thus corresponds to 
an enhanced speed down the gradient of G. However, by 
solving Eq ^fATT) with Eqs. (|A"3j) and (|A"4|) . it turns that 
when the H term is included, it accelerates the front; 
see Table |TTJ Note also that the coefficient D t is reduced 
and can even become negative in the unstable part of the 
front where C — (or at least C < C /2); see Eqs. (lA"5j) 
and (|A8|) . In that case our expression for turbulent diffu- 
sion becomes invalid and one has to include higher order 
derivatives that would guarantee stability at small length 
scales. 
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